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ABSTRACT 

Cold fronts (CFs) are found in most galaxy clusters, as well as in some galaxies 
and groups of galaxies. We propose that some CFs are relics of collisions between 
trailing shocks. Such a collision typically results in a spherical, factor « 1.4 — 2.7 
density/temperature discontinuity. These CFs may be found as far out as the virial 
shock, unlike what is expected in other CF formation models. 

As a demonstration of this effect, we use one dimensional simulations where halo 
reverberations involving periodic collisions between the virial shock and outgoing sec- 
ondary shocks exist. These collisions yield a distinctive, concentric geometric sequence 
of CFs which trace the expansion of the virial shock. 
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1 INTRODUCTION 

0\ 

Recent X-ray observations of galaxy clusters reveal various 
CN| ' phenomena in the gaseous haloes of clusters, such as merg- 
, ers, cavities, shocks and cold fronts (CFs). CFs are thought 
I to be contact discontinuities, where the density and tem- 
perature jump, while the pressure remains continuous (up 
. to projection/resolution effects). They are common in clus- 
ILJ ' ters (Markcvitch & Vikhlinin 200'S) and vary in morphology 
. !^ \ (they are often arcs but some radial or filamentary CFs ex- 
■ ist), contrast (the density jump is scattered around a value 
$^ of ~ 2), quiescence (some are in highly disturbed merging 
, regions, and some in smooth, quiet locations), and orien- 
tation. CFs have been postulated to originate from cold 
material stripped during mergers (|Markevitch et al.l [200ol . 
for example) or sloshin g of the intergalactic medium (I GM; 
iMarkevitch et all [2OO1I : lAscasibar fc MarkevitchI I2OO6I ). In 
some cases, a metallicity gradient is observed across the 
CF, indicating that it results from stripped gas, or radial 
motions of gas (see IMarkevitch fc Vikhlininll2007l . and ref- 
erences therein). Shear is often found along CFs in relaxed 
clusters, implying nearly sonic bulk fiow beneath the CF 
l|Keshet et alj|2009l ). Whether CFs are present at very large 
radii (> 500 kpc) is currently unknown because of observa- 
tional limitations. 

In what follows, we propose that some CFs are produced 
by collisions between shocks. When two trailing shocks col- 
lide, a CF is always expected to form. Its parameters can 
be calculated by solving the shock conditions and corre- 
sponding Riemann problem. Most of the scenarios in which 
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shocks are produced (quasars, AGN jets, mergers) predict 
that shocks will be created at, or near, a cluster center, 
and expand outwards (albeit not necessarily isotropically) . 
When one shock trails another, it always propagates faster 
(supersonically with respect to the subsonic downstream 
fiow of the leading shock), so collisions between outgoing 
shocks in a cluster are inevitable if they are generated within 
a sufficiently short time interval. 

In § [2] we solve for the parameters of CFs that are 
caused by a collision between two arbitrary planar shocks. 
As a demonstration of the general mechanism, we show in 
§[21 using ID spherical cluster simulation that collisions be- 
tween the virial shock and a secondary shock are expected 
to produce a distinct pattern of CFs. In §[4] we compare the 
features of observed CFs with results from our analysis, and 
list some model-specific predictions. In § [5] we summarize 
and conclude. 



2 COLLISION BETWEEN TWO TRAILING 
SHOCKS 

In this section we examine two planar shocks moving in 
the same direction and calculate the contact discontinuity 
that forms when the second shock overtakes the leading one. 
This problem is fully characterized (up to normalization) by 
the Mach numbers of the first and second shocks. Mo and 
Ml. At the instant of collision, a discontinuity in velocity, 
density and pressure de velops, corresponding to a Riemann 
problem (second case in iLandau fc Lifshitzlll959l . §93). The 
discontinuity evolves into a (stronger) shock propagating in 
the initial direction and a reflected rarefaction wave, sepa- 
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rated by a contact discontinuity which we shall refer to as 
a shock induced CF (SICF). The density is higher (and the 
entropy lower) on the SICF side closer to the origin of the 
shocks. This yields a Rayleigh- Taylor stable configuration 
if the shocks are expanding outwards. Now, we derive the 
discontinuity parameters. 

2.1 The discontinuity contrast 

We consider an ideal gas with an adiabatic constant 7. Be- 
fore the shocks collide, denote the unshocked region as zone 
0, the region between the two shocks as zone 1, and the dou- 
bly shocked region as zone 2. After the collision, regions 
and 2 remain intact, but zone 1 vanishes and is replaced by 
two regions, zone 3o (adjacent to zone 0; the outer region 
for outgoing spherical shocks) and zone 3i (adjacent to 2; in- 
ner), separated by the SICF. We refer to the plasma density, 
pressure, velocity and speed of sound respectively as p, p, u 
and c. Velocities of shocks at the boundaries between zones 
are denoted by Vi, with i the upstream zone. We rescale all 
parameters by the unshocked parameters po and po. 
The velocity of the leading shock is 
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for each zone i. Without loss of generality, we measure ve- 
locities with respect to zone 0, implying uo — 0. 

The state of zone 1 is related to zone by the Rankine- 
Hugoniot conditions, 
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pi = Po^ — f— r- — ' (3) 

(4) 

po- -. (5) 

U\ ~ vo 

The state of the doubly shocked region (zone 2) is related 
to zone 1 by reapplying equations [l]- [5] with the subscripts 
0, 1 replaced respectively by 1,2, and Mo replaced by M\. 

Across the contact discontinuity that forms as the 
shocks collide (the SICF), the pressure and velocity are con- 
tinuous but the density, temperature and entropy are not; 
the CF contrast is defined as g = pzij pzo- Regions and 3o 
are related by the Rankine-Hugoniot jump conditions across 
the newly formed shock. 
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The adiabatic rarefaction from pressure pi down to pa = 
P3i = P3o is determined by 
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The system can be solved by noting that the sum of eqs. (O 
and ((8| equals u-^—uo, fixing pzo as all other parameters are 
known from eqs. ([TjG]). 
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Figure 1. Contrast g of a CF generated by a collision between 
two trailing shocks with Mach numbers Mo and Mi, for 7 = 5/3 
(labeled contours). The color map shows the final Mach number 
A/j normalized by Mo M\. 



Finally, the Mach number of the new shock and the 
rarefacted density are given by 
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Figure [T] shows the discontinuity contrast q for vari- 
ous values of Mo and Mi, for 7 = 5/3. Typically q ^ 2, 
ranging from 1.45 for two Mach 2 shocks, for example, to 
qmax- The maximal contrast qmax depends only on the adi- 
abatic index; for 7 = 5/3 we find qmax = 2.653, achieved 
for Mo ^ 1 and M\ ~ 6.65. The possible contrast range 
of SICF s is in good agreement with the CF contrasts ob - 
served (|Markevitch fc VikhhninI l2007l : lOwers et all 120091 '). 
The figure colorscale shows the dimensionless factor / de- 
fined by Mf = fMoMi, ranging between 0.75 and 1 for 
1 < {Mo, Ml} < 100. 



2.2 Stability of Shock Induced Cold Fronts 

The stability of CFs limits the time duration over which 
they are detectable, and so is important when comparing 
CF formation models with observations. Various processes 
can cause a gradual breakup or smearing of the CF. They 
act in CFs regardless of their formation mechanism. 

Thermal conduction and diffusion of particles across 
the CF smears the discontinuity on a timescale that is set 
by the thermal velocity and the mean free path of the 
protons. The Spitzer m.f.p. A in an unmagnetized plasma 
with typical cluster densities is a few kpc, and depends on 
the thermal conditions on bot h sides of the discontinuity 
(|Markevitch fc VikhhninI |2007| ). Taking A ~ 10 kpc and a 
thermal velocity v ~ 1000 kmsec"^, and assuming that a 
CF is visible if it is sharper than Lobs ~ 10 kpc, we get a 
characteristic timescale for CF dissipation. 
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Figure 2. Fractional decline (contours) in initial CF contrast q 
induced by a collision with a shock of Mach number m arriving 
from the dense CF side, for 7 = 5/3. The vertical dashed line 
corresponds to the maximal SICF contrast, qmax- 



This result indicates that in order for shock induced CFs to 
be observable, either (i) they are formed frequently (e.g. , 
by a series of A GN bursts; see ICiotti fc Ostrikerl l2007l : 
ICiotti et al.ll2009l ): or (ii) magnetic fields reduce the m.f.p. 
co nsiderably, as s ome evidence suggests (see the discussion 
in lLazariaiill2006l '). 

Note that the heat flux driven buoyancy instability 
(HBI) tends to preferentially align magnetic flelds perpen- 
dicular to the heat flow in regions where the temperature 
decreases in the direction of gravity. This effect could re- 
duce radial diffusion within t he inward cool i ng regions in 
the c ores of cool core clusters (|Quataertll2008l : iParrish et al.l 
Across a CF transition there is a sharp temperature 
drop towards the inner side, in the direction of gravity. Heat 
flow could induce HBI on the scale width of the CF, aligning 
the magnetic fields parallel to the discontinuity and dimin- 
ishing radial diffusion. This possibility needs to be addressed 
further by detailed numerical magneto-hydrodynamic simu- 
lations. 

CFs could be degraded by subsequent shocks that sweep 
outwards across the CF. However, the CF contrast is only 
slightly diminished by this effect, by < 4% for g < 2 and 
a shock Mach number m < 2, as shown in fig. [21 Such 
passing shocks seed Richtmyer-Meshkov instabilities which 
could cause the CF to break down. Although less efficient 
than Rayleigh- Taylor instabilities, these instabilities oper- 
ate regardless of the alignment with gravity. The outcome 
depends on m and on any initial small perturbations in the 
CF surface. However, the instability is suppressed if the CF 
becomes sufficiently smoothed. 

It is worth noting that KHI, which could break down 
CFs formed through ram pressure stripping and sloshing. 



does not play an important role in SICFs because there is 
little or no shear velocity across them. On the other hand, 
the stabilizing alignment of magnetic fields caused by such 
shear (jMarkevitch fc Vikhlinin. ,2007i ) is also not expected 
here. 



3 A CLUSTER SCALE REVERBERATION 
MODE 

At the edge of galaxy clusters there are virial shocks with 
typical Mach numbers ~ 30 — 100, heating the gas from 
~ 10* K to ~ 10^ K by converting kinetic energy into ther- 
mal energy. The rate of expansion of a virial shock is set on 
average by the mass fl ux and velocity of the infalling ma- 
terial (see for example iBertschingej [l985l ) . We flnd in our 
simulations that while this expansion is steady on average, 
the shock sometimes oscillates, moving faster or slower than 
the steady state expansion. This creates periods dominated 
by alternating halo compression and expansion. During peri- 
ods of enhanced compression in the outer halo regions, when 
the shock slows down, compression is sent inwards and is re- 
flected through the center, gradually steepening into an out- 
going shock. When this shock collides with the virial shock, 
a situation corresponding to that described in § [21 bouncing 
the virial shock into another cycle of reverberation, creat- 
ing an SICF and sending inwards a steep rarefactioiij. This 
rarefaction marks the end of the compression period of the 
next cycle setting the oscillation period to approximately 
the sound crossing time, in and out across the cluster. The 
combined effect is a long-lived, coherent "breathing" mode. 
Such "breathing" oscillations a re observed in the 3D galactic 
and cluster halo simulations of iKeres fc Hernguistl (|2009l 'FI. 

Using ID spherical simulations of th e formation of clus- 
ters from initial cosmic perturbations (|Birnboim fc Dekej 
we test the formation and evolution of associated 
SICFs. The simulation includes baryonic shells, dark mat- 
ter shells, and angular momentum support. The oscillatory 
mode in the simulation has not been intentionally excited, 
and results from stochastic inner core (~ 0.05Rvir) vibra- 
tions as external dark matter shells with low angular mo- 
mentum interact gravitationally with the core. Physically, 
any waves, shocks or gravitational perturbations (mergers or 
AGNs) will contribute to the mode, as well as non-smooth 
accretion rates that would cause the virial shock to vibrate. 
A more detailed examination of this is beyond the scope of 
this letter (note, however, that the strength of the SICF de- 
pends only weakly on the parameters of the secondary shock 
and is always close to 2.) 

Fig. [3l shows the evolution of a 3 X 10" M© halo un- 
til redshift z — (an adiabatic version of those described 
in Birnboim & Dekel 2009, m prep.). The expanding virial 
shock is traced by a jump in the velocity of the Lagrangian 
shells (red, thin lines), and by the shock-flnding algorithm 
that is based on entropy increase along Lagrangian shells 
(blue region). The outgoing, secondary shocks that reach the 
virial shock periodically are also clearly visible. The virial 
shock bounces outwards each time it is hit by a secondary. A 

^ A movie of the evolution of radial profiles in time is available 
at http : //www. cf a. harvard. edu/~ybirnboi/SICF/sicf .html 
^ Dusan Keres, (private communication). 
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Figure 3. Evolution of a galaxy cluster from cosmological ini- 
tial perturbations. The final mass of the cluster at z = is 
3 X IO^^Mq. The radius and time of Lagrangian shells (every 
25*'' shell) are plotted in red thin lines. Shocks are traced by 
large Lagrangian derivatives of the entropy {dlnS/dt > .1 Gyr~^, 
blue dots), and CPs are traced by their large entropy gradients 
(91n5/<91nr > 0.5, green dots). Rarefaction waves can be seen as 
small motions of the Lagrangian shells (illustrated by the dashed 
magenta curve, manually added based on a time series analysis). 
Their trajectory is approximately the reflection of the preceding 
outgoing compression, time inverted about the last secondary- 
virial shock collision. 



rarefaction wave is reflected inwards (illustrated by a dashed 
curve), and a spherical SICF is left behind. The SICFs are 
visible through the entropy gradients (green region), and, in 
the absence of diffusion and instabilities, persist until z — 0. 
As expected, the interaction of the discontinuity with sub- 
sequent shocks ( i|2.2l fig. [2]) does not considerably reduce 
the CF contrast. Fig. [4] shows the entropy, pressure, den- 
sity and temperature profiles of this simulation at z = 0. A 
series of SICFs (marked with blue vertical lines) is visible, 
seen as entropy jumps with continuous pressure. The den- 
sity jumps by roughly a factor of q ~ 2, as expected, with 
the temperature dropping accordingly. 

The SICFs formed in this simulation are almost static, 
and roughly logarithmically spaced. The calculations pre- 
sented in Figures [3]|4] are adiabatic. When cooling is turned 
on in the absence of feedback, this simulation suffers from 
overcooling, creating an overmassive BCG of 2x10^2Mq and 
star formation rates of 100 Moyr"^, and the luminosity ex- 
ceeds the Lx—T relation. The excitation of the basic mode, 
and periodic CFs, are observed in all these simulations. The 
simulation presented here was performed with 2, 000 bary- 
onic and 10, 000 dark matter shells. When the resolution is 
lowered to 250 baryonic shells, the reverberation amplitude 
and periods are essentially unchanged, indicating that the 
results are well converged. 




> 



Radius^pc] 

Figure 4. Thermodynamic profile of the simulated cluster shown 
in Fig. Q ai z = 0. Top: entropy (red, solid, left axis) and pres- 
sure (green, dashed, right axis). Bottom: density (red, solid, left 
axis) and temperature (green, dashed, right axis). CFs are marked 
with blue dashed vertical lines. 



OBSERVABILITY OF SHOCK INDUCED 
COLD FRONTS 



lOwers et al.l (|2009l ) present high quality Chandra observa- 
tions of 3 relaxed CFs. They all seem concentric with re- 
spect to the clust er center, and spherical in appearance. 
lOwers et al.l (|2009l ) interpret these as evidence for sloshing, 
and find spiral characteristics in 2 of them. We argue here 
that some CFs of this type could originate from shock col- 
lisions. In the SICF model, CFs are spherical, unlike the 
truncated CFs formed by ram pressure stripping or slosh- 
ing, so no special projection orientation is required. On the 
other hand, SICFs do not involve metallicity discontinuities, 
observed in some of these CFs. The observed contrasts of 
these CFs are quite uniform, with best fit values in the range 
2.0 — 2.1 (with ~ 20% uncertainty) in all three, as expected 
for SICFCl. 

The SICF formation model is entirely distinct from ram 
pressure stripping and sloshing CFs, and it makes many dif- 
ferent predictions. 

Morphology. SICFs are quasi-spherical around the 
source of the shocks. A galactic merger defines an orbit 
plane, and the corresponding perturbation (stripped mate- 
rial or center displacement) will cr eate CFs parallel to this 
plane. In some sloshing scenarios (jAscasibar fc MarkevitchI 
120061) the CFs extend considerably above the plane, but 
would never appear as a full circle on the sky; an observed 
CF "ring" would strongly point towards an SICF. The sta- 
tistical properties of a large CF sample could thus be used 
to distinguish between the different scenarios. 

Amplitude. SICFs have distinct entropy and density 
contrasts that depend weakly on shock parameters; q is typ- 



^ Note that collisions with outgoing shocks could diminish the 
contrast towards q ^ 2, regardless of CF origin. A larger, more 
complete sample of CFs is needed to determine the origin/s of 
the relaxed population. 
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ically larger than ~ 1.4 (assuming Mi ^ 2) and is always 
smaller than qmax = 2.65. 

Extent. If the cluster reverberation mode is excited, 
SICF radii should be approximately logarithmically spaced. 
Any shock that expands and collides with the virial shock 
will create SICF at the location of the virial shock, far be- 
yond the core. The external SICFs are younger, and would 
appear sharper owing to less diffusion and collisions with 
secondary shocks. Deep observations capable of detecting 
CFs at r ~ 1 Mpc are predicted to find SICFs (fig-U). Such 
distant CFs occur naturally only in the SICF model. 

Plasma diagnostics. Shocks are known to modify 
plasma properties in a non-linear manner, for example 
by accelerating particles to high energies and amplify- 
ing/generating magnetic fields. The plasmas on each side 
of an SICF may thus differ, being processed either by two 
shocks or by one, stronger shock. This may allow indirect de- 
tection of the CF, in particular if the two shocks were strong 
before the collision. For example, enhanced magnetic fields 
below the CF may be observable as excess synchrotron emis- 
sion from radio relics that extend across the CF, in nearby 
clusters, using future high-resolution radio telescopes (LO- 
FAR, SKA). 



and rep r oduce some of the CF features discussed in 
lOwers et al.l l|2009l '). In particular, we suggest that CFs in 
relaxed clusters, with no evidence of mergers, shear, or chem- 
ical discontinuities, may have formed by shock collisions. We 
list the properties of SICFs that could distinguish them from 
CFs formed by other mechanisms. The SICF model predicts 
quasi-spherical CFs which are concentric about the cluster 
center, with contrast g ~ 2, and possibly extending as far 
out as the virial shock. An observed closed (circular/oval) 
CF could only be an SICF. In the specific case of cluster 
reverberation, a distinct spacing pattern between CFs is ex- 
pected. It may be possible to detect them indirectly, for ex- 
ample as discontinuities superimposed on peripheral radio 
emission. 

Shocks originating from the cluster center naturally oc- 
cur in feedback models that are invoked to solve the over- 
cooling problem. They are also formed by mergers of sub- 
structures with the BCG. Thus, SICFs should be a natural 
phenomenon in clusters. Further work is needed to assess 
how common SICFs are with respect to other types of CFs, 
and to characterize inner SICFs that could result, for exam- 
ple, from collisions between offset AGN shocks. The proper- 
ties of reverberation in 3D will be pursued in future work. 



5 SUMMARY AND DISCUSSION 

Our study consists of two parts. The first is a general discus- 
sion about CFs that form as a result of a collision between 
two trailing shocks. The second part describes a specific re- 
verberation mode of galaxy and cluster haloes that has been 
identified in ID simulations, and also seen in 3D. One as- 
pect of these reverberations is periodic collisions between 
the virial shock and outgoing secondary shocks, resulting in 
SICFs - a special case of the mechanism discussed in the 
first part of our paper. 

We have shown that when shocks move in the same 
direction they collide and generate a CF. The density con- 
trast across the CF is calculated as a function of the Mach 
numbers of the two shocks. It is typically larger than 1.4 (if 
Af > 2), and is always smaller than qmax = 2.65. The dis- 
continuity is smeared over time by diffusion, at a rate that 
depends on the unknown nature and amplitude of magnetic 
fields. CFs are susceptible to heat-fiux-driven buoyancy in- 
stability (HBI), which could align the magnetic field tangent 
to the CF and potentially moderate further diffusion. SICFs, 
like all other CFs, are subject to Richtmyer-Meshkov insta- 
bilities from subsequent shocks passing through the cluster. 
Such collisions reduce the CF contrast until it reaches g ~ 2. 
Unlike most other CF models, an SICF is not expected to 
suffer from KHI. 

In § [S] using a ID spherical hydrodynamic code to 
evolve cosmological perturbations, we demonstrate that a 
reverberation mode exists in the haloes of galaxies and clus- 
ters and causes periodic collisions between the virial shocks 
and secondary shocks that produce SICF. The simulated 
SICF contrast is consistent with the theoretical predictions 
of § m A more thorough investigation of this potentilly im- 
portant mode, including its stability in 3D is left for future 
works. We use it here as a specific scenario for SICF forma- 
tion. 

The predicted properties of SICFs are presented in 



ACKNOWLEDGEMENTS 

We thank M. Markevitch for useful discussions. YB ac- 
knowledges the support of an ITC fellowship from the Har- 
vard College Observatory. UK acknowledges support by 
NASA through Einstein Postdoctoral Fellowship grant num- 
ber PF8-90059 awarded by the Chandra X-ray Center, which 
is operated by the Smithsonian Astrophysical Observatory 
for NASA under contract NAS8-03060. 



REFERENCES 

Ascasibar Y., Markevitch M., 2006, ApJ, 650, 102 
Bertschinger E., 1985, ApJS, 58, 39 
Birnboim Y., Dekel A., 2003, MNRAS, 345, 349 
Ciotti L., Ostriker J. P., 2007, ApJ, 665, 1038 
Ciotti L., Ostriker J. P., Proga D., 2009, ApJ, 699, 89 
Keres D., Hernquist L., 2009, ApJ, 700, LI 
Keshet U., Markevitch M., Birnboim Y., Loeb A., 2009, 
ArXiv e-prints 

Landau L. D., Lifshitz E. M., 1959, Fluide Mechanics. Perg- 

amon press, 1959 
Lazarian A., 2006, ApJ, 645, L25 

Markevitch M., Ponman T. J., et al. 2000, ApJ, 541, 542 
Markevitch M., Vikhlinin A., 2007, Phys. Rep., 443, 1 
Markevitch M., Vikhhnin A., Mazzotta P., 2001, ApJ, 562, 
L153 

Owers M. S., Nulsen P. E. J., Couch W. J., Markevitch M., 

2009, ApJ, 704, 1349 
Parrish I. J., Quataert E., Sharma P., 2009, ApJ, 703, 96 
Quataert E., 2008, ApJ, 673, 758 



